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Experiments on several novel superconducting compounds have observed oscillations of the specific 
heat when an applied magnetic field is rotated with respect to the crystal axes. The results are 
commonly interpreted as arising from the nodes of an unconventional order parameter, but the 
identifications of nodal directions are sometimes controversial. Here we show with a semiclassical 
model calculation that when the magnetic field points in the direction of the nodes, either minima 
or maxima can occur in the specific heat depending on the the temperature T and the magnetic 
field H. An inversion of the angular oscillations takes place with respect to those predicted earlier 
at low temperature by the nodal approximation. This result, together with the argument that the 
inversion takes place based on an approximation valid at moderate fields, indicates that the inversion 
of specific heat oscillations is an intrinsic feature of nodal superconductors. 

PACS numbers: 



I. INTRODUCTION 



Initial information on the symmetry of the order pa- 
rameter in newly discovered superconductors is often pro- 
vided by power laws in the temperature dependence of 
thermodynamic and transport properties 1 . However, the 
exact nodal structure cannot be determined by these 
techniques because the same exponents may correspond 
to order parameters with different symmetries. Phase- 
sensitive experiments like the tricrystal experiments per- 
formed on cuprates 2 are the most definitive, but are tech- 
nically challenging and require high quality samples. A 
simpler technique, proposed in Ref. 3, which is not phase 
sensitive, but provides information on the distribution of 
nodes on the Fermi surface, is a measurement of the spe- 
cific heat in the presence of a magnetic field, H, which is 
rotated with respect to the crystal axes of the sample. It 
was shown in that work that, if the effect of the supcrflow 
due to vortices on the quasiparticle spectrum was treated 
semiclassically, via the Doppler energy shift 4,5 , the low- 
temperature specific heat of a superconductor with line 
nodes acquires an oscillatory dependence on the field ori- 
entation. In particular, the low-temperature specific heat 
coefficient 7 ~ T\/~H was shown to have minima when 
the magnetic field points in the direction of the nodes 
in the order parameter, and maxima for the field along 
the antinodes, where the gap is maximal. Ref. 3 and its 
extension, Ref. 6 made use of two additional approxima- 
tions. First, they replaced the Doppler energy shift of 
a quasiparticle with momentum k near one of the nodes 
by its value at the node, k = k„. This so-called nodal 
approximation had been shown to work well compared 
to the full semiclassical evaluation in Ref. 5, providing 
T,Eh *C Ao- Here Ao is the maximum gap over the 
Fermi surface, and Eh ~ Aq\/ H/H C 2, where H C 2 is the 
upper critical field, is a magnetic energy scale (see Sec. 
II). Second, making predictions for low temperatures and 
fields, T,Eh <C Ao, the authors of Refs. 3,6 used a form 



of the order parameter linearized in the vicinity of the 
nodal points. 

By now there have been several experimental tests of 
these ideas 7-13 , and the observations have been generally 
consistent with theoretical expectations. At the same 
time, the assignment of the nodal directions in several 
materials remains controversial. In CeCoIns the mea- 
surements of the anisotropy of the specific heat 9 and the 
thermal conductivity 1 ' 15 appear to give contradictory 
results for the gap structure. In Sr2Ru04 the specific 
heat oscillations were observed to invert as the temper- 
ature and field was varied, i.e. the minima and maxima 
as a function of angle changed places 10 ' 13 . 

Such an inversion was never found in earlier theoret- 
ical calculations within the semiclassical approach, even 
though it was found in theoretical work employing other 
techniques, see below. This clearly poses a problem: if 
the technique is to be useful as a way to "map" out the 
nodal structure, it is necessary to be able to predict if 
and when such inversions will occur; otherwise the nodes 
in the order parameter may be assigned to incorrect lo- 
cations in momentum space. 

Numerical solution of the Bogoliubov-dcGennes equa- 
tions yielded an inversion in the anisotropy of the density 
of states (DOS), N(u, 9), between the field applied in the 
nodal and antinodal directions 16 , but this was argued to 
reduce, rather than invert, the specific heat anisotropy. 
Recently Vorontsov and Vekhter 17,18 considered the limit 
H c i <C H < H C 2 by extending the method of Brandt, 
Pesch and Tewordt 19,20 , and found an agreement with 
semiclassical method at low T, H, but an inversion of the 
specific heat oscillations over a large part of the T — H 
phase diagram. Since the approximation they used was 
tailored for moderate to high fields, their determination 
of location (or, indeed, the existence) of the inversion line 
could be questioned. 

In this paper we demonstrate that this inversion is a 
generic feature of the specific heat in unconventional su- 
perconductors. We revisit the semiclassical approach to 
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the in-plane specific heat oscillations of a quasi-two di- 
mensional d-wave superconductor, but relax the approxi- 
mations which led in Rcfs. 3,6,21 to simple analytical forms 
for C(T; 9) at low temperatures. We are thus forced to do 
a numerically demanding evaluation of the entropy and 
specific heat, which involves a 2D k-summation and a 2D 
averaging over a vortex lattice unit cell. Our model is not 
restricted to low temperatures, and is valid (within limits 
discussed below) in the low field regime H c \ < H <C H C 2- 
It may thus be considered to complement the results of 
Refs. 17,18 . We find that the anisotropy of the specific 
heat is inverted as the temperature is increased. This 
result demonstrates that the inversion phenomenon is 
robust across the phase diagram for unconventional su- 
perconductors, and provides an important caveat to the 
interpretation of the stationary points in specific heat 
oscillations. 



II. THE SEMICLASSICAL APPROACH 

In superconductors with line nodes, the vortex core 
contribution to the low-energy density of states is smaller 
than that from the quasiparticles outside the core 
region 4 . In systems with short coherence length, such 
as cuprates and heavy fermion materials, the dominance 
of the extended quasiparticle states is even more pro- 
nounced since the spacing of the energy levels in the vor- 
tex core is large, and only few such levels (if any exist at 
all) are occupied at low temperature. Bulk quasiparticles 
in the vortex state are excited from the pair condensate 
moving with the superflow around each vortex; hence the 
effect of applied magnetic field can be simply described by 
Dopplcr-shifting the spectrum of extended quasiparticle 
states according to the local value of the superfluid veloc- 
ity, v s (r). This approximation is valid at H <C H C 2, when 
the vortices are far apart and v s (r) varies slowly on the 
scale of the superconducting coherence length. Using this 
approach, Volovik predicted 4 that the density of states 
at the Fermi level in the vortex state of a superconductor 
with line nodes varies as N(lj = 0; H) oc s/H. Specific 
heat measurements on high-T c cuprates verified the \[H 
field dependence 22-24 , which played an important role 
in the identification of d-wave symmetry. It should be 
noted, however, that the result holds for any gap with 
line nodes, including, for example, some p-wave states. 
The \f~H dependence is modified by the presence of dis- 
order to HlogH as shown by Kiibert and Hirschfeld 25 , 
but remains nonlinear and qualitatively similar to the 
pure case. 

The semiclassical approximation provides a conceptu- 
ally transparent and tractable approach to the effect of 
the vortex lattice on low-lying quasiparticle states. Its 
validity has been questioned by several authors, who at- 
tempted a more accurate description of the effect of both 
applied field and supcrcurrents on the quasiparticle spec- 
trum. Franz and Tcsanovic 26 introduced a singular gauge 
transformation that takes into account both the super- 



current distribution and the magnetic field, and mapped 
the full quantum-mechanical problem onto that of nodal 
Dirac fermions interacting with effective scalar and vec- 
tor potentials that are periodic in the unit cell of the vor- 
tex lattice. They discovered significant differences from 
the quasiclassical theory in the quasiparticle excitations 
at very low energy, as is to be expected since the semi- 
classical approach works only for large quantum num- 
bers. We do not discuss the details of this breakdown 
here, since it was done by Knapp et al. 27,28 , who explic- 
itly compared the quantum-mechanical and semiclassical 
results for the density of states, and found that small 
differences between the two approaches begin to appear 
below a crossover scale which is exponentially small in 
the Dirac cone anisotropy vf/v2- Here vf is the Fermi 
velocity, and V2 = dA/dku is the gap "velocity" at the 
node (fcy is the momentum component along the Fermi 
surface). Since in real materials this ratio is large, and 
in real samples the presence of impurity scattering and 
the disorder in the vortex lattice smear out the energy 
structure on small scales, we assume that for purposes 
of comparison with the measurements, the semiclassical 
description is adequate. 

Damn et al. 29 investigated in detail the comparison 
of the simple single vortex Doppler shift approach with 
the solution of the quasiclassical microscopic Eilenbcrgcr 
equations using several techniques including the Brandt, 
Pesch, and Tewordt approximation 19,20 . At distances 
from the vortex of the order of the coherence length, 
the Doppler shift method is quantitatively inadequate be- 
cause of core state contributions or scattering resonances, 
but these effects have little qualitative impact on the cal- 
culation of thermodynamic properties. We therefore pro- 
ceed with the simplest Doppler-shift analysis for a single 
vortex unit cell, in order to make the qualitative point 
which is the main result of this work. 

For concreteness, we consider a quasi-two dimen- 
sional superconductor with d-wave symmetry. This sit- 
uation is most closely realized in cuprates, possibly the 
heavy fermion 115 compounds 14,30 , and potentially also 
the newly-discovered oxypnictide materials 31 . We con- 
sider the field applied parallel to the conducting plane, 
H || ab, and assume that the system is sufficiently three- 
dimensional, so that at the fields of interest the structure 
of the mixed state resembles the Abrikosov vortex lat- 
tice penetrating a stack of weakly coupled 2D planes 32 . 
In this case the vortex superflow field v s is three dimen- 
sional, different on different planes within the vortex unit 
cell. For weak interlayer coupling we consider a circular 
in-plane Fermi surface, and therefore account only for 
quasiparticles with momenta k || ab. Deviations from 
Fermi surface isotropy and effects of multiple Fermi sur- 
face sheets have been studied e.g. in Refs. 33 and 34 , but 
do not affect the qualitative conclusions wc wish to draw 
here. 

The single-particle Green's function in the presence of 
a superflow velocity field v s is obtained by Doppler shift- 
ing the quasiparticle states with energy uj and momentum 
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k 5 (we use units with ft = &g = 1). 



G(k, r,u n ) 
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(iu) n - v s (r) • k) 2 - £ 



k 



(1) 



where Lu n is the fermionic Matsubara frequency, £k is 
the band energy measured with respect to the Fermi 
level, Ti are Pauli matrices in particle-hole space, and 
Ag = Ao cos 20, with <j) the azimuthal angle on the Fermi 
surface. 

We consider the magnetic field, H, applied in the 
ab plane, and approximate the supcrflow by that of an 
anisotropic 3D superconductor in the London model. 
The contours of constant v s are then elliptical due to 
anisotropy of the penetration depth, X a (, ^ X c . After 
rescaling the c-axis by z' = zX a b/X c , the superflow is 
cylindrically symmetric and the Doppler shift for quasi- 
particles at the Fermi surface is given by 3 ' 6 
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where ip is the winding angle of the superfluid velocity 
in real space, 9 is the angle between H and the a-axis, 
and r is the rescaled distance from the center of the vor- 
tex. Note that we explicitly exclude from consideration 
quasi- ID Fermi surfaces such as those shown by Tanaka 
et al. 34 to lead to C(H; 9) oscillation inversion at low 
temperatures. 

From Eqs. (1) the local density of states is 
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Here No is the normal state density of states. The net 
DOS per volume is found by spatially averaging N(u>, r) 
over a unit cell of the vortex lattice containing one flux 
quantum, $o- After rescaling the c-axis the unit cell 
area increases by a factor X c /X a b, and therefore flux 
quantization dictates that the quasiparticles now expe- 
rience the effective field H* = (X a b/X C )H. In the new 
coordinates the radius and the area of such a cell are 
Rh = sj &o/irH* and Ah = nR H respectively. Intro- 
ducing polar coordinates, r = Rh (pcos ip, psimp), we 
find for the average field-dependent density of states 



N(u,H,T) = - f P dp f , 

T JO JO 



di/jN(uj,r) . 
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The density of states depends on the angle between 
the field and the crystal axes. At low energies the domi- 
nant contribution to the local DOS is from the near-nodal 
regions, <f> ~ <f> n = ir/A ± nn/2, with n = . . . 3, and os- 
cillations in N(u>, H, T) appear since the Doppler shift at 
a given node vanishes when the field is aligned with that 
node. In the nodal approximation with linearized order 
parameter, the analytical form of the density of states 



was obtained 3 ' 5 ' 6 , see next Section; however, as pointed 
out before, this approximation does not give the reversal 
of the anisotropy (see below) . 

In all of the following we carry out the full summation 
over the momenta and the real space coordinates. We 
use the density of states to compute the entropy, 



S = -2 / dojN(uj)[(l - /(w)) log(l - /M) 

J — oo 

+/(a/)log/H] > (5) 

where f(uj) = (exp(cj/T) + is the Fermi function. 
To obtain the the heat capacity at constant volume we 
differentiate, 
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At low T and H, when the gap varies weakly with tem- 
perature, the temperature derivative acts only on the 
Fermi function in Eq.(5), and the specific heat is given 
by the simple form 
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Here we consider the full range of temperatures below 
the transition, and hence evaluate the specific heat from 
Eq. (6). We approximate the temperature dependence of 
the order parameter, Ak(T), according to the BCS weak- 
coupling ansatz appropriate to a circular Fermi surface 35 



A(T) = A cos(20) tanh 
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A V 3 7C(3) v T 
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where Ao = 2.14T C is the gap maximum. 

In the absence of impurity scattering, there are two im- 
portant low energy energy scales in the problem: the tem- 
perature T and the magnetic energy, or typical Doppler 
shift Eh = vp/Rn- To satisfy the requirements of the 
semiclassical approach, we consider only Eh «C Ao, but 
temperature to vary over the entire range T < T c . 



III. DENSITY OF STATES IN PLANAR FIELD 

In the nodal approximation of Ref. 3, the quasiparticlc 
momentum k is replaced by its value k n at each of the 4 
nodes, which are then summed over. The residual density 
of states at the Fermi level, Eq. (4) then becomes 
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where the angular variation is given by P{9) = 
max (| sin 0|, I cos #|). Eq. (7) then yields the linear in 
temperature specific heat at low T with the slope 
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At low energies, Eh,oj <C Ao, the density of states in 
the nodal approximation with linearlized order parame- 
ter takes the form 3,6 



N(u,H,T) = 1 
N 2 
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where £i = E H \ sin(7r/4 — 6)\, E 2 = E H \ cos(7r/4 
and the scaling function F is given by 5 



F(y) 



y [1 + l/(2j/ 2 )] , 
(1 + 2y 2 ) arcsiny 



if y > i; 

%\/i -y 2 /y*, if y < i- 

(12) 

Note that in the limit y — > oo, ^(y) — * y, such that 
N(u),H, T) in (4) recovers the isotropic low-w d-wave 
density of states lu/Aq. Thus in the nodal approxima- 
tion with linearized order parameter the specific heat os- 
cillations are washed out at higher temperatures, but the 
method cannot generate specific heat oscillation inver- 
sions. 
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FIG. 1: (Color online) Comparison of the nodal approxima- 
tion Eq.(ll) for the density of states N(u,H) versus the full 
density of states at T = 0.001T C and E H = 0.2T C . solid 
line: nodal approximation for H along antinode. Dashed line: 
nodal approximation for H along node. Red circles: full eval- 
uation for H along node. Black triangles: full evaluation for 
H along antinode. The insert magnifies the low-frequency 
range. 

In Fig. 1, we now compare the density of states of a 
d-w&ve superconductor obtained from a complete evalu- 
ation of Eq. (4) with the nodal approximation with lin- 
earized order parameter Eq. 11. In agreement with all 
previous work the density of states for field in nodal and 
antinodal directions is strongly anisotropic for lo < Eh, 
and the anisotropy is washed out at higher energies. 
However, the difference between JV(w,H, T) for the two 
directions of the field reappears at energies of order Ao , 
the energy scale absent in the versions of the calculation 
with linearized gap. 

For the field along the antinode i?j|k a „ (6 = 0,7r/2...) 
the density of states continues to be sharply peaked at 
Ao, as in the absence of the field. This peak is largely 



due to the quasiparticlcs moving along the field, which do 
not experience the Doppler shift and see the full maximal 
gap. In contrast, for the field along the node, if ||k„ 
(6 = 7r/4, 37r/4, . . . ), the DOS has broad features around 
Ao ± Eh- As a result, the density of states for the field 
along the nodal direction begins to exceed that for field 
along the antinode at the lower shoulder, and the DOS 
anisotropy is inverted, relative to that at low energies. 

We emphasize that the absence of the structure near 
the gap edge in N(u),H, T) in previous work is a conse- 
quence of gap linearization, and not the nodal approxi- 
mation. We observe features similar to those depicted in 
Fig. 1 when replacing k — ► k„ in the Doppler shift, but 
keeping the full variation of the order parameter around 
the Fermi surface. Of course, the nodal approximation is 
not accurate at the energies of order Ao, and hence we 
continue with the full evaluation of the DOS. 

The anisotropy in the density of states that we com- 
puted differs from that found in Refs. 17,18 , where no- 
ticeable inversion of the anisotropy in N(ui, H, T) oc- 
curred already at relatively low energy, as a result of the 
competition between Doppler shift and vortex scattering. 
Remarkably, we find that even within the semiclassical 
method, when the lifetime remains infinite and the scat- 
tering on the vortices in neglected, the anisotropy in the 
density of state is still reversed, albeit at higher energies 
of order of the gap maximum. 

The possibility of an inversion of the specific heat 
anisotropy is clear from Fig. 1. In Eq. (7) the density 
of states is convoluted with the temperature-dependent 
weighting function, peaked around u> ~ 2.5T. At a given 
field, as the temperature is increased, more weight in the 
kernel shifts to higher energies, where the anisotropy in 
the DOS is opposite to that at low co. Whether an in- 
version then occurs at higher T depends primarily on 
whether the kernel has sufficient weight in the expo- 
nential tail at energies to ~ Ao — Eh- Since the DOS 
anisotropy is small in magnitude above this crossing en- 
ergy, the specific heat must be calculated numerically. 

In Figure 2 we plot the density of states as a field 
sweeps through the ab-plane. We clearly see that for 
low frequency the node and minimum in the oscillations 
coincide, and that at higher frequencies this is no longer 
the case. The eventual inversion of this pattern results in 
the angular density of states having a maximum at the 
gap node for higher frequencies. 



IV. SPECIFIC HEAT OSCILLATIONS 

Numerical differentiation of the entropy is computa- 
tionally intensive due to the high accuracy required in 
finding S in Eq. (5). To illustrate the precision of our 
calculations in Fig. 3, we show the numerically evaluated 
specific heat at H = and at Eh = 0.4T C for the field 
along the a-axis over a wide temperature range both be- 
low and above T c . The normal state specific heat above 
T c is C N = j N T = 2A tt 2 T/3, where is the Fermi 
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FIG. 2: DOS N(9, uj) vs. magnetic field angle 9 at E h = 0.2T C 
and T = 0.05T C for equally spaced frequencies, 0.025Ao apart, 
from uj = to 0.175Ao. The solid line is the result of the nodal 
approximation for the uj = DOS from Eq. (9). 



level DOS per spin. 




FIG. 3: Temperature dependence of the specific heat (at con- 
stant volume) normalized to the specific heat in the normal 
state Cv/Cn, where Cat = JnT. We show both the exact 
(symbols), and the asymptotic low-T (dotted line) behavior 
in the absence of the field, and compare it to the numerically 
determined Cv/Cn at Eh = 0.4T C (dashed line). There is no 
reduction of T c with field in our approach. 

To test our numerical evaluations we first check the nu- 
merical results against the asymptotic low-temperature 
specific heat in zero field. For our model of a 2D e?-wave 
superconductor with a circular Fermi surface at T <C T c 
we find 



Cy 
NnT 



7o 



T 



7o = a- 



dx 



A J cosh 2 (z/2) 



7Tl8C(3) 



(13) 
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which agrees with the numerically determined slope. The 
Volovik effect is manifested in Fig. 3 in the finite offset 



of C/T (Eq. (10) in the presence of magnetic field Ejj 
0AT c . 
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FIG. 4: Upper panel: C(T, H)/C(T, 9 = 0) vs. 9 for a set of 
equally spaced temperatures, every 0.02T C , from to 0.18T C , 
Eh = 0.2T C . Lower panel: same for a set of equally spaced 
temperatures from to 0.08 T c every 0.02T C , E h = 0.05T C . 

We proceed to evaluate the specific heat for fixed Eh 
and several temperatures as a function of the field angle 
9. Fig. 4 shows that the inversion of the DOS anisotropy 
found in Sec. Ill indeed are sufficient to lead to the in- 
version of the specific heat oscillations a characteristic 
temperature, T inv . In Fig. 4 the fourfold oscillations are 
clearly visible at low T, and minima occur for H along 
nodal directions as anticipated. However, at higher tem- 
peratures, an inversion in the pattern of oscillations is 
evident. In Fig. 5 we examine the anisotropy by plotting 
the difference between the specific heat for the field along 
the nodal and the antinodal directions, which identifies 
the temperature at which inversion occurs. 

Since in our approach the inversion is due to the sen- 
sitivity of the specific heat to the changes in the DOS 
within energy range ~ Eh of the gap edge, increasing 
both the magnetic field and the temperature initially en- 
hances the amplitude of the inverted oscillations. In- 
creasing the field brings the anisotropy inversion down 
in energy in Fig. 1; raising the temperature increases the 
contribution of the high energy regions in C/T. While 
the amplitude of the inverted (relative to those at T = 0) 
oscillations in Cy is small, it is of the same order of mag- 
nitude as that observed experimentally 7, 9-12 , and found 
theoretically for a quasi-two-dimensional system at mod- 
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FIG. 5: Amplitude of specific heat oscillation de- 
fined by difference of Cv between node and antinode, 
[CV(0)-CV(f)]/CV(0) for E H = 0.2T C (solid line) and 
E H = 0.05T C (dashed line). 



crate fields 17,18 . Of course, when H (T) approaches the 
upper critical field (the transition temperature) respec- 
tively, the gap in the spectrum closes, and the oscillations 
vanish; we cannot, however, reliably comment on the evo- 
lution of the anisotropy in this regime within the semi- 
classical method. At the same time it follows from our 
analysis that the amplitude of the inverted oscillations 
has a maximum at intermediate fields and temperatures, 
also in agreement with Refs. 17 ' 18 . Therefore our results 
connect well with those obtained by a different technique. 



V. CONCLUSIONS 

In this paper we have calculated the specific heat of a 
two-dimensional d— wave superconductor in an external 
magnetic field using the semiclassical treatment of the 
effect of the vortex lattice on the quasiparticle spectrum. 
In contrast to previous work utilizing the nodal approx- 
imation with linearized order parameter, we carried out 
a full numerical evaluation of the density of states and 
the entropy for a wide range of fields and temperatures. 
Our main finding is that the sign of the oscillations of the 
specific heat as a function of the field orientation, i.e. the 
difference between CV for the field along a nodal direc- 
tion and along the antinode, depends on the temperature 
and field strength. We confirmed that at low tempera- 
tures and fields the specific heat has a minimum when 
the field is along a nodal direction. However, as H and T 
are increased, minima of the specific heat begin to occur 
for the field along the gap maxima, i.e. an inversion of 
the oscillation pattern occurs. Note that while we consid- 
ered a system with well-defined quasiparticle states, it is 
reasonable to believe that scattering due to impurities or 
vortex lattice disorder will merely smear the anisotropy 
on both sides of the inversion line. 

Our calculations provide a bridge connecting the semi- 
classical theory at low fields H « H C 2 with the re- 
sults of the extended Brandt-Pcsch-Tcwordt approxima- 



tion 17 ' 18 , where the inversion was first found. The latter 
approach is in principle valid at H < H C 2, but has been 
shown to provide remarkably good agreement with semi- 
classical predictions down to low fields, up to log correc- 
tions 18 . To compare the results of the two approaches 
explicitly, we extend our calculation to higher fields, ac- 
count for the field dependence of the gap amplitude via 
A(T, H) = A(T) V /1 - H/H c2 , and determine the inver- 
sion line. This provides a direct comparison with the 
results of Ref. 18, where the same field dependence was 
assumed for a circular Fermi surface identical to that 
considered above. In Fig. 6 we plot the approximate 
crossover scales for both the Brandt-Pesch-Tewordt and 
semiclassical theories. There is a remarkably good quali- 
tative agreement for the behavior of the inversion line up 
to moderate fields, where the ranges of validity of the two 
approaches may reasonably be assumed to overlap. This 
establishes a phase diagram for when the specific heat is 
expected to have minima or maxima at the gap nodes. 
Taken together, these results strongly suggest that the 
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FIG. 6: Phase diagram H/H C 2 vs. T/T c for specific heat oscil- 
lations of a d-wave superconductor in a magnetic field. The 
shaded regions indicate that the specific heat C(T, H) has 
a minimum when the field H points in the nodal direction, 
whereas the white regions indicate inverted oscillations, where 
minima correspond to field along antinode. The gray shaded 
regions are results from Ref. 18 using the Brandt-Pesch- 
Tewordt framework, whereas the blue regions represent the 
region of the H — T plane where minima correspond to nodes 
within the semiclassical theory. The phase diagram compar- 
ing both approaches was determined using the assumption 
H/H c2 = (Eh/A ) 2 , valid up to a prefactor of order unity. 
Insert: blowup of the low T and H region. Note the field scale 
is given in terms of .Eff/An for easier comparison with other 
results in this work. 

inversion of the specific heat oscillations with tempera- 
ture is a general feature for all nodal superconductors, 
and establish the approximate location of the inversion 
line. 

Consequently, identification of the nodes in the gap 
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via the specific heat measurements is not as straightfor- 
ward as the original semiclassical results suggested, and 
depends on where the experiment is done in the field- 
temperature plane. While the inversion of the oscillations 
in related experiments on the anisotropy of the heat con- 
ductivity of nodal superconductors has been well estab- 
lished 15, 36-38 , an inversion in the heat capacity measure- 
ments has not yet been searched for systematically. Ex- 
perimental identification of the inversion line in the T-H 
plane would be an important step towards further estab- 
lishing the heat capacity measurements as the method of 
choice for determining the nodal directions in the bulk, 
and our paper provides theoretical foundation for such a 



search. 
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